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19 Key Points 

20 • Interhemispheric gradient can serve as a proxy to constrain CCI 4 emissions. 

21 • A minimum 30 Gg/yr CCI 4 emissions is necessary to reconcile the observations. 

22 • The likely lifetime for CCI4 is 25-36 years, longer than the current estimate. 
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Abstract 


Carbon tetrachloride (CCI4) is a major anthropogenic ozone-depleting substance and greenhouse 
gas and has been regulated under the Montreal Protocol. However, atmospheric observations 
show a very slow decline in CCI4 concentrations, inconsistent with the nearly zero emissions 
estimate based on the UNEP reported production and feedstock usage in recent years. It is now 
apparent that there are either unidentified industrial leakages, an unknown production source of 
CCI4, or large legacy emissions from CCI4 contaminated sites. In this paper we use a global 
chemistry climate model to assess the budget mystery of atmospheric CCI4. We explore various 
factors that affect the global trend and the gradient between the Northern and Southern 
hemispheres or interhemispheric gradient (IHG): emissions, emission hemispheric partitioning, 
and lifetime variations. We find a present-day emission of 30-50 Gg/yr and a total lifetime ~25- 
36 years are necessary to reconcile both the observed CCI4 global trend and IHG. 

Index terms 

Constituent sources and sinks; Troposphere: constituent transport and chemistry; Air/sea 
constituent fluxes. 


Page 3/25 


42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 


1. Introduction 


Carbon tetrachloride (CCI 4 ) is primarily used as a feedstock or processing agent for 
chlorinated species, but has been used extensively as a cleaning agent and as a solvent in the past 
[CTOC Report, UNEP , 2011]. CCI 4 is recognized as both an ozone-depleting substance (ODS) 
and a greenhouse gas. As of 2008, CCI 4 accounted for about 11% of total tropospheric chlorine 
[ WMO , 2011]. The ozone depletion potential (with respect to CFC-1 1) is 0.82 \WMO, 2011] and 
it has a 100-yr global warming potential of 1,400 [WMO, 2011]. In 1987, Article 2 of the 
Montreal Protocol (MP) included regulations of CCI 4 under Annex B Group 2. CCI 4 production 
and consumption were eliminated for developed countries in 1996 under the amendments to the 
MP. Developing countries (i.e., Article 5 countries) were allowed some phase down production 
and consumption until fully banned in 2010. CCI 4 continues to be legally used as a contained 
feedstock, e.g. for the production of hydro fluorocarbons (HFCs), since feedstock uses are not 
regulated by the Montreal Protocol. 

The primary sinks for CCI 4 include photolysis loss in the stratosphere, degradation in the 
ocean and the soil [SPARC, 2013]. The current best estimate of total lifetime (x) for CCI 4 is 25 
years [SPARC, 2013], relatively unchanged from the WMO [2011] assessment (26 years). A best 
estimate of the atmospheric partial lifetime (T atm 0 s) for CCI 4 is 44 years [SPARC, 2013]. The best 
estimates of partial lifetimes due to the ocean sink (x 0C ean) and the soil sink (x so ii) are 
approximately 81 (71-167) years and 195 (108-907) years [SPARC, 2013], respectively. 

The MP controls have led to declining CCI 4 levels in our atmosphere at a rate slightly greater 
than 1% per year [WMO, 2011]. Under the MP Article 7, each Party provides CCI 4 data to the 
Ozone Secretariat on production, imports, exports, feedstock amounts, and amounts destroyed. 
The current CCI 4 bottom-up emissions estimate from the MP parties based on reported 
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production and feedstock usage (bottom-up estimate) was zero after 2007 [ WMO , 2011]. There 
are also no known substantial stocks of CCU in existing equipment or storage containers; thus, a 
x ~ 25 years would imply an annual decrease of 4% per year rather than the observed 1%. The 
atmospheric CCI4 observations and the current total lifetime can be used to derive a top-down 
estimate of global emissions. This top-down emission estimates for 2007 and the following years 
were upward of 50 Gg per yr (Gg/yr) [WMO, 2011]. This very large difference of emission 
estimates is equivalent to approximately ~ 1,600 railroad tank cars of liquid CCI4 lost each year. 

The difference between the top-down and bottom-up emission estimates suggests that there is 
an unreported source of CCI4. Recent work by Fraser et al. [2013] suggests that emissions from 
contaminated soils, toxic waste treatment facilities, and possibly chloro-alkali plants could be 
contributing 10-30 Gg/yr. De Bias et al. [2011] also observed excess CCI4 above the background 
in Bilbao, Spain (similar to the Fraser et al. observations near Melbourne, Australia), and 
attributed this to an unidentified source near the measurement site. Odabasi [2012] found that 
mixing of bleach with surfactants or soap could form CCI4, but global emissions from this source 
have not been estimated. Emissions from CCI4 feedstock uses are highly uncertain [TEAP, 2011], 
but have been estimated to be approximately 0.5% of the total feedstock used (equivalent to 5 
Gg/yr for 2011 production) [ Miller and Batchelor, 2012]. None of these potential sources alone 
can fully explain the 50 Gg/yr discrepancy between the top-down and bottom-up emission 
estimates. 

In this paper, we use available source and sink data in global and box models to test the 
compatibility of the existing emission and lifetime estimates with CCI4 mixing ratio observations. 
In particular, we present the most likely emission and lifetime scenarios that close the current 
gap and best reconcile the observed trend in concentrations and their difference between the 
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Northern hemisphere (NH) and Southern hemisphere (SH) (i.e., the inter-hemispheric gradient or 
IHG). 

2. Models and Simulations 

Models. We pair a 3-Dimensional (3-D) Chemistry Climate Model (CCM) with 1-box and 2- 
box models to examine the atmospheric budget of CCU. The global 1-box model used in this 
study is the same model used in recent WMO Ozone Assessments, and is described in detail in 
Velders and Daniel [2013]. The global 1-box model has been used to derive top-down emission 
estimates for long-lived ODSs using best-estimate lifetimes and observed surface mixing ratios 
as constraints. The 2-box hemispheric model is developed from the global 1-box model and 
simulates both the long-term global trend and Northern- Southern hemispheric differences in 
atmospheric concentrations of CCI4. We assume ocean and soil losses in two hemispheres are 
scaled exactly with ocean and soil area. The IHG is generated by the asymmetry in hemispheric 
emission fraction and the ocean and land surface fractions in the two hemispheres. The 3-D 
CCM used here is the NASA GEOS Chemistry Climate Model ( GEOSCCM) Version 2, which 
couples the GEOS-5 GCM [ Reinecker et al., 2008] with a detailed stratospheric chemistry 
module [ Douglass and Kawa, 1999]. A comprehensive evaluation of several CCMs over the 
1960-2005 period shows that the GEOSCCM agrees well with observations for many of the 
meteorological, transport-related, and chemical diagnostics [Eyring et al, 2006]. Of particular 
relevance to this study, GEOSCCM represents well the mean atmospheric circulation as 
demonstrated by its realistic age-of-air, and further, realistic loss and atmospheric lifetimes for 
long-lived ODSs [Waugh et al., 2007; Douglass et al., 2008; Chipperfield et al., 2014]. The 
model also features realistic inter-hemispheric transport and reproduces well the observed IHG in 
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previous flux-based simulations of major long-lived ODSs [Liang et al., 2008; Chipperfield et al., 
2014]. A detailed description of the box models and the GEOSCCM is in Appendix A. 

While the 3-D GEOSCCM is used to understand how various processes impact CCI4 
concentrations in a more realistic modeled atmosphere, the box models provide top-down 
emissions and lifetime estimates that are consistent with the observed surface mixing ratios. 
Pairing the simple box models with the 3-D GEOSCCM greatly enhances the computational 
efficiency of choosing emissions and corresponding lifetimes for the 3-D model simulations. The 
box models also provide useful conceptual tools to examine the sensitivities of CCI4 mixing ratio 
observations to sources, sinks, and the atmospheric inter-hemispheric distributions of these 
quantities (2-box hemispheric model). 

3-D Simulations. The global emissions in the 3-D simulations are top-down emission 
estimates consistent with the observed atmospheric CCI4 decline when using the global 1-box 
model. CCI4 is run with flux boundary conditions, using geographically resolved surface 
emissions originally described in Xiao et al. [2010]. Five simulations are performed (Table 1). 
The baseline run, Run A, is a 53-yr simulation with the SPARC 2013 photochemistry, soil and 
ocean lifetime recommendations, and the corresponding top-down emission derived using the 
global 1-box model. To better understand the factors that influence the CCI4 budget, we conduct 
four additional model simulations (1995-2012), Runs B-E, with varying lifetimes, global 
emission and emission distributions. Runs B-E are initialized with the January 1995 Run A 
initial conditions. 

3. Results 

3.1 Discrepancy between bottom-up and top-down emissions estimates 
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We use the long-term surface observations of CCI4 made by the National Oceanic and 
Atmospheric Administration - Global Monitoring Division (NOAA-GMD) [ Montzka et a/., 

1999; Thompson et al., 2004] to derive the top-down emission estimates in the global 1-box 
model. The GMD dataset is a combination of in situ and flask samples, all based on Gas 
Chromatography - Electron Capter Detector (GC-ECD) analysis [Hall et ah, 2011]. The global 
mean atmospheric CCI4 is decreasing at a mean rate — 1.1 ppt per year since 1995 (Figure 1A). 
With a x - 25 years, this suggests a slow decrease in emission from ~ 80 Gg/yr from the late 
1990s’ to ~ 55 Gg/yr in the early 2010s’ (Figure 1C). Reported industrial production and 
feedstock usage of CCI4 suggests a rather sharp decrease in CCI4 emissions from 100 Gg/yr in 
1999 to near-zero emissions after 2007 (Figure 1C) if all CCI4 is emitted in the year in which it 
was produced, as is generally expected. Such a drastic emission reduction rate is inconsistent 
with the observed CCI4 decline. 

We use the deviation of CCI4 surface mixing ratios from the linear decay line to estimate 
year-to-year changes in annual emissions. After removing the least-squares linear fit (2000-2012) 
from each NOAA GMD station, we apply a 25-month ^-amplitude Gaussian low-pass filter to 
the observations (Figure IB). The filtering reveals 3 periods of change: 1) from 1995-2005 a 
continuous increase in CCI4 anomalies across all stations (mean — +0.2 ppt/yr) with an increase of 
0.6 ppt from 2003 to 2005 (+0.3 ppt/yr), 2) from 2007-2011 a decrease of about 1 ppt (-0.25 
ppt/yr), and 3) an anomalous jump of ~0.6 ppt beginning in about 2012 (+0.6 ppt/yr). Based on 
the global 1-box model estimate, a 1 ppt change in atmospheric mean CCI4 is equivalent to - 25 
Gg/yr change in emissions. Thus, these observed anomalies imply: 1) - +8 Gg/yr anomaly 
emissions in period 1 from 2003-2005, 2) a — 6 Gg/yr extra emission decrease between 2005- 
2008, and 3) an anomalous increase of + 15 Gg/yr in 2012. These year-to-year changes in 
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observed CCI 4 anomalies are inconsistent with the bottom-up emissions estimate from reported 
production and consumption (Figure 1C). 

These inconsistencies between the mean trend and year-to-year emission variations and the 
observed CCI 4 changes suggest that the bottom-up emissions estimate from reported production 
and consumption are likely incorrect. 

3.2 The inter-hemispheric gradient 

It has long been recognized that IHG is a qualitative indicator of emissions for long-lived 
chemical compounds [Lovelock et a/., 1973]. The results from our 3-D model simulations show a 
compact linear correlation between the model annual IHG and the annual global emissions for all 
individual runs (R = 0.92-0.98 for Runs A-D) (Figure 2a). Of course, changing the hemispheric 
emission ratio affects this ratio, as do changes to the distributions of loss between the 
hemispheres (see discussion below). The collective correlation coefficient between the two 
variables from all four runs is 0.96, despite the various emissions and lifetimes used in each 
model run. This implies that IHG can be used as an empirical proxy to quantitatively infer global 
emissions. 

It is important to mention that model results suggest that the global mean IHG calculated 
using all grid points in each hemisphere is different than that calculated using only model values 
at the NOAA GMD stations, due to biased sampling using only station data (Appendix B). The 
model global mean IHG is ~1.2 ppt higher than the IHG calculated using the model grid point 
values sampled at the GMD stations. This difference remains rather constant between 1995-2012. 
As a result, we apply a +1.2 ppt correction factor to the IHG calculated using the NOAA GMD 
station observations (1.5 ± 0.3 ppt). This corrected IHG is referred to as the NOAA GMD- 
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inferred IHG (IHGgmd) in the rest of the paper, which ranges between 2. 3-3.0 between 1995- 
2012 . 

3.3 3-D model simulations: Budget constraints from trend and inter-hemispheric gradient 

We use the global trend and the IHG as two independent constraints in the 3-D GEOSCCM. 
All model runs, except Run E, are designed to reproduce the observed CCI4 global trend. The 
model IHG is then compared with IHGgmd to assess various emission and lifetime scenarios that 
best reconcile the observed IHG. 

There are many factors that contribute to the CCI4 IHG, including global emissions, 
hemisphere emission fractions (EFhemis), and soil and ocean loss rates. Despite a large range of 
emission strengths, total and partial lifetimes used, runs A-D yield very similar IHG-emission 
regression slopes (0.049-0.058 ppt/Gg yr' 1 ) (Figure 2a). These regression lines also show similar 
zero-emission intercept points at 0.59-0.64 ppt (Figure 2a), the likely IHG that can be explained 
by ocean and soil losses alone. 

Of the range of parameters explored in the 3-D simulations, global emission strength plays 
the dominant role in determining the IHG. Baseline Run A employed the highest emissions and 
yielded the highest IHG between 3-5 ppt for 1995-2012, ~ 50% higher than IHGgmd- This 
suggests that the mean ~ 64 Gg/yr emissions estimate in Run A is likely biased high. To test this 
emission level, Run B was employed with the lowest mean emission considered of ~ 35 Gg/yr. 
Run B had the smallest IHG between 2. 3-3.0 ppt for 1995-2012, agreeing well with IHGgmd- 
This is not surprising as Run B was designed from the 2-box model with global emissions that 
would reproduce the observed gradient, albeit a corresponding lifetime increase to ~ 36.5 years 
was necessary to match the long-term trend. The two runs with intermediate emissions (mean ~ 
50 Gg/yr) produce intermediate IHGs (2. 6-3. 7 ppt for Run C and 2. 8-4.0 ppt for Run D). 
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In addition to changes in global emissions, Run C is designed to explore the sensitivity of 
IHG to changes in EFhemis- We change EFhemis from the baseline 94%NH:6%SH used in Runs A 
and B to 88%NH: 12%SH. Run C shows a slightly smaller IHG/emissions regression slope 
(0.049 ppt/Gg vs. 0.053 ppt/Gg in Run A). This implies, to match the observed gradient, higher 
global emissions are needed if one assumes a larger fraction of emission resides in the SH. The 
result of this magnitude of repartitioning emissions into the SH is a relatively small reduction of 
the IHG, in comparison to the global emission strength. 

The oceanic loss also affects the IHG. Run D features a latitude dependent ocean loss with a 
faster degradation in the SH. Faster SH ocean loss rates lead to a slight increase in the IHG (=< 
10%). Again, this ocean loss rate impact on the IHG is small, in comparison to the global 
emission strength. 

Run E is a special run in which we used global emissions consistent with the observed IHG, 
but lifetimes were kept the same as in SPARC [2013]. The Run E CCU decreased at ~ 2.2 ppt/yr, 
double the observed rate. This suggests that while an average global emission ~ 35 Gg/yr is in 
better accordance with the observed IHG and the bottom-up estimate, there is a large discrepancy 
between this estimate and the current best estimate x of 25 years for closing the global CCI4 
budget. 

3.4 What impacts the inter-hemispheric gradient? - Insights from the 2-box hemispheric 
model 

The 2-box hemispheric model yields a similar a strong linear relationship between the IHG 
and global emission as the 3-D GEOSCCM. With the inter-hemispheric exchange timescale 
(tinterhemis) set at ~1.7 years, the regression slope of IHG vs. emissions from the 2-box 
hemispheric model reference calculation (0.05 ppt/Gg, Figure 3 red symbols) agrees well with 
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the 3-D model hemispheric mean IHG from runs using the same EFhemis and partial lifetimes. 

The zero-emission intercept from the 2-box model is, however, only ~0.1 ppt, much lower than 
the ~0.6 ppt value from the 3-D model. This is likely due to missing atmospheric processes that 
could impact the IHG, e.g. stratosphere-troposphere exchange differences between the NH and 
SH, or a simplified constant Tinterhemis- Nevertheless, the similar IHG/emissions relationship 
between the 3-D global model and the 2-box hemispheric model makes it possible to employ the 
2 -box hemispheric model to explore the important factors that detennines the IHG. 
Mathematically, IHG can be approximated using the following equation: 

IHG = a + b x Emission (1) 

Where a (unit of ppt) is the zero-emission intercept point and b (unit of ppt/Gg yr' 1 ) is the 
regression slope on the scatter diagram. The value of a is a function of x OC ean, x SO ii, EFhemis, and 
likely STE as well in the 3-D model, and b is mainly decided by EFhemis- 

We vary EFhemis, T 0C ean, T S0 ii, and Tatmos to illustrate how changes in each variable regulate the 
IHG (Figure 3). All 2-box hemispheric model calculations shown here use consistent emissions 
and lifetimes that reproduce the observed global CCU trend between 1995-2012. While the top- 
down emission estimates are fixed for any specified lifetime input, the IHG can vary 
significantly with changing EFhemis (Figure 3a). A 100% NH emission (b = 0.066 ppt/Gg yr' 1 ) 
requires - 30-40 Gg/yr global CCI4 emissions to reproduce the IHGgmd- Decreasing the NH 
release fraction greatly increases the total emissions necessary to capture the IHGgmd (e.g, 
b=0.026 ppt/Gg yr' 1 for 70% release fraction in the NH, which requires emissions >75-115 
Gg/yr). This 70% scenario is unlikely as the needed emissions greatly exceed the bottom-up 
estimate - approaching the peak emissions - 120 Gg/yr in the 1970s’ and 1980s’ before CCI4 was 
regulated by the MP. 
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Changing partial lifetimes also impacts the IHG. Increases of T 0C ean (he., less loss, see Fig. 

3b) decreases the IHG by decreasing the hemispheric contrast in ocean surface losses, reflected 
in the value of a in Eq. (1). As x oce an increases from 80 years to 240 years, a decreases from 0.21 
ppt to 0.05 ppt. Increases of x so ii (Fig. 3c) increase the IHG, and a increases from 0.08 ppt to 0.22 
ppt as Tsoii increases from 200 years to 1000 years. Changing T atm0 s has little impact on a (Fig. 3d). 
Overall, these impacts are small. However, increases in x OC ean, T SO ii, and T atm os can affect the 
calculated IHG to a greater extent through increases in total lifetime and the implied necessary 
decreases in emissions to match the observed trend. Hence, the 2-box hemispheric model 
calculations, similar to the 3-D model results, indicate that global emissions and EFhemis play a 
dominant role in controlling the IHG with partial lifetimes associated with ocean, soil, and 
atmospheric losses contributing only minor modifications. 

3.5 The likely emissions and lifetime scenarios 

The 2-box hemispheric model, even with the least number of necessary processes considered 
in the present form, contains more unknown variables (x, global emissions, EFhemis) than 
constraints (observed trend and IHG). Hence, it is inadequate to uniquely close the CCfi budget 
mystery. Alternatively, we use the least root mean square deviation (RMSD) approach with our 
current best understanding of emissions and loss processes to rule out the unlikely EFhemis 
scenarios and, furthennore, to infer an optimal total lifetime and emission scenario for each 
likely EFhemis, as illustrated in Figure 4. 

The RMSD of the calculated IHG suggests that the NH fractional emission is unlikely to be 
less than 80%, as a 70% or 76% NH emission fraction would yield an optimal total lifetime of ~ 
20 years. This not only implies a much larger gap in total emissions needed to reconcile with the 
current bottom-up emissions estimate, but also significant decreases in partial lifetimes, both of 
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which are difficult to accommodate. The RMSD approach suggests that the likely NH emission 
fraction is between 80%-100%. This range implies average emission ~ 30-50 Gg/yr for 1995- 
2012 and optimal x ~ 25-36 years, equal or longer than the current best estimates. A smaller NH 
fractional emission release requires a corresponding scenario of more global emissions and 
shorter x to match the observed gradient, e.g. emission ~ 50 Gg/yr and x ~ 25 years for a 80% 

NH emission fraction. 

Assuming 100% emissions reside in the NH (b=0.066 ppt/Gg yr' 1 ) and using a mean IHGgmd 
~ 2.7 ppt and a ~ 0.6 ppt from the 3-D model in Eq. (1), we are able to determine the minimum 
mean global emissions necessary to match the atmospheric CCI 4 observations is ~ 32 Gg/yr 
between 1995-2012 and ~ 29 Gg/yr for the 2010s’. However, this requires a lifetime of 36 years, 
much longer than the 25 years current estimate. In the 3-D model Run B, we tested this 36-year 
lifetime by increasing x a tmo S to ~ 62 years and x OC ean to ~ 160 years and leaving x so ii unchanged. An 
increase in Run B’s x atm0 s from 47 years to 62 years requires a ~ 60% reduction in the CCI 4 
photolysis rate which greatly exceeds the lab-measured 15-20% cross section uncertainty range 
\Rontu Carton et al., 2010; SPARC, 2013]. When comparing Run B’s CCI 4 with two limited 
balloon profiles, this leads to a model high-bias in the critical stratospheric photolysis loss region 
(10-70hPa) (Figure A3). Keeping x atm 0 s unchanged, ax ~ 36 years means the summed lifetime 
against ocean and soil losses is > 150 years, which requires increases of both x oce an and x so ii to the 
upper limit of the current best estimates (167 years for x ocean and 907 years for x so ii). The above 
points to a need of re-evaluation of the current best estimates of partial lifetimes to address the 
possibility of a longer x. 

4. Summary and discussions 
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CCI4 was increasing in the atmosphere until the early 1990s and is now in decline [WMO, 
2011]. This decline is a result of the regulations imposed by the Montreal Protocol. The decline 
during the 1990-2006 period was caused by a decrease in emissions and removal from the 
atmosphere via loss processes [WMO, 2011; SPARC, 2013]. 

The current CCI4 downward trend should be primarily determined by the lifetime, because 
bottom-up emissions after 2007 are estimated to be near zero. However, from 2007-2012, the e- 
folding time-scale of the decrease in the CCI4 observations is about 66 years. This slow 
downward trend cannot be reconciled with our current best estimate of the ~ 25 years (implied 
top-down emissions ~ 55 Gg/yr between 2007-2012) derived from comprehensive chemistry- 
climate models and atmospheric, ocean, and soil observations [SPARC, 2013]. 

Comparisons of mixing ratio observations with year-to-year variations of bottom-up 
emissions estimates cannot be understood either. Much larger year-to-year CCI4 fluctuations are 
expected from the current time-series of bottom-up emissions estimate than are actually observed 
in the atmospheric abundances. For example, between 1996 and 1998 the estimate of CCI4 
emissions jumped by about 80 Gg. Cumulatively, this 3-year period would have resulted in a 12 
ppt increase in CCI4 concentrations with respect to the slow decreases due to atmospheric, soil, 
and ocean losses. Observations filtered to show these shorter time scale changes reveal 
approximately a 1-2 ppt change of CCI4. Such discrepancies indicate potential serious flaws with 
CCI4 emission reports. 

Using a fully-coupled chemistry-climate model with a state-of-the-art photochemical loss 
scheme for CCI4, along with current estimates of CCI4 oceanic and soil sinks, we have performed 
a series of model simulations to constrain the atmospheric budget of CCI4. The inter-hemispheric 
gradient (IHG) of CCI4 has been qualitatively used in the past to in fer emissions of long-lived 
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ozone-depleting substances. Our COM results suggest that the IHG provide useful information 
for quantitatively estimating CCI4 emissions. We exploited the global trend and IHG as two 
independent constraints in this study to evaluate possible explanations for the discrepancy 
between top-down and bottom-up emission estimates. 

The near-zero emissions from the UNEP reported production and feedstock usage in the 
recent years cannot be reconciled with these model estimates. At a minimum, a present-day 
global emission of 30 Gg/yr is required to reproduce the observed CCI 4 trend and IHG. It is most 
likely that ~80%-100% of the total global emissions are released in the NH. The likely range of 
global emissions and total lifetime (x) associated with the above range of NH emission fraction 
remains large, varying from global emission ~ 30 Gg/yr and x ~ 36 years for a 100% NH 
emission fraction to global emission ~ 50 Gg/yr and x ~ 25 years for a 80% NH emission fraction. 
In the majority cases, this implied x needed to reconcile the observed trend and IHG, is longer 
than the current best estimate lifetime (25 years). This would necessitate longer atmospheric, 
ocean and/or soil partial lifetimes than the current best estimates. Our findings point to the need 
of a more accurate bottom-up emissions estimate and/or lifetime estimate to close the CCI 4 
budget mystery. Alternatively, information on fractional emission estimate from the two 
hemispheres, which is not well quantified currently, can also help in narrowing the likely range 
of global emission and lifetime. 
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407 Table 1 . A description of the five 3-D GEOSCCM CCfi simulations used in this work. 



Description 

Simulation 

Period 

Average 

emission 

1995-2012 

(Gg/yr) 

Hemispheric 

emission 

fraction 

NH:SH 

Lifetime 

X 

(yr) 

Partial lifetimes 

(yr) 

^atmos 

Xocean 

T-soil 

Run A 

Baseline simulation 

1960-2012 

64 a 

94%:6% 

25.8 

47 

80 

200 

Run B 

Decreased ocean loss, 
Decreased atmospheric loss 
forced by reducing the 
photolysis rate 

1995-2012 

35^ 

94%:6% 

36.5 c 

62 c 

160 c 

200 

Run C 

Repartitioning of emissions 
into the SH with reduced global 
emissions. 

1995-2012 

50 e 

88%: 12% 

30.7 d 

47 

160 d 

200 

Run D 

1 As Run C, but with latitude- 
dependent ocean loss rates with 
faster degradation in the 
Southern Hemisphere. 

1995-2012 

50 

88%: 12% 

29.5 

47 

135 

200 

Run E 

Same lifetimes as in Run A and 
same emissions as in Run B. 
This simulation does not match 
the observed CCfi decline. 

1995-2012 

35 

94%:6% 

25.8 

47 

80 

200 


408 a The global 1-box model top-down emissions estimate for x ~ 25.8 years. 


409 b IHG-based annual emissions calculated using the average of IHGgmd (section 3.2) and the IHG 

410 from the Advanced Global Atmospheric Gases Experiment (AGAGE) network [ Prinn et al., 

411 2000] (Appendix B). 

412 c x, Xatmos, and x oce an are determined using the global 1-box model in the forward mode with the 

413 IHG-based emissions and the observed global trend. 

414 d For Run C, x oce an is determined using the 2-box hemispheric model by matching the IHGgmd- 

415 e The global 1-box model top-down emissions estimate for x ~ 30.7 years. 

416 1 The latitude-dependent ocean loss rates used are 1/288 yr' 1 for 45-90°N, 1/222 yr' 1 for 0-45°N, 

417 1/122 yr" 1 for 0-45°S, 1/75 yr" 1 for 45-90°S. The relative strength of latitude-dependent loss rates 

418 are provided by Shari Yvon-Lewis (personal communication) and then scaled to give an ocean 

419 partial lifetime of 135 yrs. 
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Figure 1. A) The observed CCI4 mixing ratios at the NOAA GMD stations (color symbols) and 


the global mean values (solid black line). B) Same as A) but for mixing ratio anomalies. Note 


that some 2012-2013 data are preliminary, and have not yet undergone final calibration. C) CCI 4 


emissions derived from atmospheric measurements (red and blue shading) and bottom-up 


potential emissions estimated from production data (green lines). The potential emissions 


estimate was derived from the difference between total CCI4 production reported to UNEP (black 


line) and feedstock (line labeled). Red shading indicates the top-down emissions estimate from 


the global 1-box model using a total lifetime range of 25-36 yr. Blue shading indicates the IHG- 


scaled emissions using IHGgmd and an empirically derived scaling factor with a range of 1/0.05- 


1/0.06 Gg yr'Vppt (section 3.5). 
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Figure 2. (a) The scatter diagram of model mean inter-hemispheric gradient vs. global annual 
emissions used in each model run between 1995-2012. Each symbol represents one annual- 
averaged value. The dashed lines (same color as the symbols) show the regression slope for each 
run. The gray shaded region indicates the IHGgmd between 1995-2012. (b) The global mean 
CCI4 mixing ratios from the NOAA GMD stations (thick black line) and model runs A-E. 
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Figure 3. The red plus symbols on all panels show the scatter diagram of annual mean IHG vs. 
global emissions between 1995-2012 from a reference calculation from the 2-box hemispheric 
model. This reference calculation is equivalent to the 3-D model Run C, with a hemispheric 
emission fraction of 88%NH:12%SH, T OC ean = 160 yrs, T S oif=200 yrs, andT atm os = 47 yrs. Each 
symbol represents one annual-averaged value. The gray shaded region indicates the range of 
IHGqmd between 1995-2012. The groups of blue squares on each panel show the 2-box 
hemispheric model sensitivity results by varying one of the input variables, (a) EFhemis with 
regression slope b (ppt/Gg yr' 1 ) shown in parenthesis, (b) T OC ean, (c) x so ii, and (d) x atm os, with 
respect to the reference calculation. The dashed lines (same blue color as the symbols) show the 
regression slope for each corresponding 2-box hemispheric model calculation. 
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Figure 4. The root-mean-square-deviation (RMSD) of 2-box hemispheric model calculated IHG 
(using the IHGgmd as references) as a function of total lifetime for each assumed hemispheric 
emission fraction (70%-100%). For each hemispheric emission fraction line, we highlight in red 
symbols the optimal total lifetime that yields the least RMSD. 
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